Field-domain spintronics in magnetic semiconductor multiple quantum wells 
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We develop a theory of non-linear growth direction transport in magnetically doped II- VI com- 
pound semiconductor multiple-quantum-well systems. We find that the formation of electric field do- 
mains can be controlled by manipulating the space dependence of the band electron spin-polarization, 
using its exchange coupling to local moments. We emphasize the importance of band electron spin 
relaxation in limiting the strength of these effects. 
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I. INTRODUCTION 



The giant niagnetoresistance effect in magnetic metal 
multilayerstJ occurs because of the coupling of external 
magnetic fields to band electron spins through their col- 
lective spin polarization. The utility of this effect for in- 
formation storage and field sensing devices has increased 
interest in exploring related spin-dependent transport 
properties in both ferromagnetic and paramagnetic semi- 
conductors. At the same time, progress in the homo 
and hetero epitaxy of magnetically doped semiconduc- 
tors is creating new possibilities for engineered material 
geometries in whichjiew spin-dependent transport effects 
are likely to occup.au Large band electron spin polar- 
izations iSan occurQ in diluted magnetic semiconductors 
(DMSs),eI for example in II- VI compounds with Mn sub- 
stituted on the group II site. Among the striking phe- 
nomena already demonstrated in Mn doped II- VI quan- 
tum structures are magnetically tunable quantum well 
barriers and interwell couplings,Q spin-dependent dynam- 
ics of polarized excitona.in spin supcrlattices,Ll optically- 
probed spin coherence,S and the injection of highly po- 
larized|-spin currents into GaAs/AlGaAs light emitting 
diodes.Ll In addition, n-doping of wide-gap II- VI mag- 
netic semiconductotjauantum wells (Zni_3;_j,Cd2;MnySe) 
has been achieved ,E3 yielding two-dimensional electron 
gas (2DEG1 systems that are exchange coupled to mag- 
netic ionsO The exchange interaction, Jsd, between s- 
electrons in the conduction band and Mn++ S — 5/2 lo- 
cal moments, results in band spin splittings larger than 
hujc, the Landau level splitting.-IThe spin-splittings can 
reach values as high as 20 meV.o In fact, complete spin 
polarization can be achieved in quantum wells at rela- 
tively low magnetic fields (~ 1 T).E3 

Although the study of electronic transport prop- 
erties in DMS heterostructure systems is still in 
its initial stages, interesting predictions have already 
been madcllJ^Ej The present study is motivated in 
part by the recent growth of a modulation-doped 
ZnSe/(Z«,Cd,Mn)Se multiple quantum weU (MQW) 
system.ta In non-magnetic MQW systems, growth di- 



rection transport phenomenology is enriched by an in- 
terplay between charge accumulation and resonant inter- 
well tunneling effects that results in the formation of elec- 
tric field domains. In this paper we report on a theory 
of the influence of exchange coupling with Mn ion spins 
on electric field domain formation and on the sensitiv- 
ity of this influence to spin-relaxation rates within the 
quantum wells. 

The formation of electric field domains is the hallmark 
of dc-biased transport in weakly-coupled semiconductor 
superlattices. Spontaneously generated inhomogeneities 
in the spatial distribution of voltage drops, were pro- 
posed as aruepcplanation for overall conductance {G{V)) 
oscillations,c_l discovered in pioneering growth direction 
MQW transport studies. This early hypothesis was latei. 
confirmed by direct photoluminiscence measurements. E3 
Subsequently, highly doped GaAs/AlGaAs superlattices 
that present sawtooth-like current-voltage (I-V) charac- 
teristics in the negative differential conductance (NDC) 
region were studied in detail. Along branches of the saw- 
tooth two approximately constant electric field regions 
develop in the sample, separated by a layer of accumu- 
lated electrons. In the following, we follow common usage 
in referring to these layers with higher 2D electron den- 
sity as monopoles. This non-equilibrium configuration 
enables resonant tunneling between ground and excited 
subbands in the high field region, minimizing the total re- 
sistance of the superlattice. Increases in external voltage 
in this regime lead to sharp decreases in current, followed 
by discrete jumps of the monopole region from a well to 
its upstream neighbor, extending the high field domain 
over an additional period and increasing the current. 

Past work has studied the dependence of domain for- 
mation and evolution on magnetic fields applied along 
the growth direction and on far-infrared radiation wave- 
length. In the former case, the formation of Landau levels 
and scattering between them_.introduces a new voltage 
scale for domain formation.Ej In the latter case, pho- 
tonic sidebands sustain the formation of the electric field 
domains.Ej The inclusion of the electronic spin in the 
study of perpendicular transport in MQWs can also be 
expected to alter electric field domain formation physics. 



and exchange coupling to Mn spins should make it pos- 
sible to tune these effects with relatively weak exter- 
nal magnetic fields. This is the possibility that we ex- 
plore at greater length in this paper. In all these cases 
the electron-electron interaction, although small in the 
ground state in comparison with typical energies of the 
system, can not be neglected since it is the Poisson equa- 
tion relating charge accumulation to field variations that 
is at the heart of field-domain formatimi. For example it 
permits the experimentally observedEil multistability of 
distinct stationary physical states at a fixed bias voltage. 
Nonetheless, the mean-field Hartree approximation is suf- 
ficient to capture this physics in typical samples. The 
non-linearity of the current versus voltage relationship 
between neighboring quantum wells, coupled with the 
non-locality of electron-electron interaction effects leads 
to transport equations that can be solved only numeri- 
cally, and also to results that are sometimes difficult to 
interpret. 

In this paper we deal with the formation of electric 
field domains in II- VI MQW systems with one or more 
(II,Mn)VI quantum wells. The main ingredients of our 
self-consistent theoretical model are: (i) a theory for the 
tunneling current between two spin-polarized 2DEGs; (ii) 
a continuity equation that accounts for relaxation of non- 
equilibrium spin populations; (iii) a relationship between 
the up and down chemical potentials and their densities; 
(iv) the application of simple Hartree mean field theory 
to account for the Coulomb interaction; and (v) a mean- 
field theory for the interaction between 2DEG electrons 
and Mn spins whose average polarization is very sensi- 
tive to external magnetic fields at low temperature. We 
shall demonstrate that new features appear in the I-V 
curve that depend on temperature and Mn spin concen- 
tration, and explain why spin bottlenecks turn out to 
have a strong influence in the instability regions. 

The paper is organized as follows. In Section I| our 
theoretical model is thoroughly explained. Sec. |lll is de- 
voted to a discussion of technical details important for 
the numerical integration of the rate equations that de- 
scribe the time dependent charge, spin, and current dis- 
tributions. In Sec. IV we give numerical results and dis- 
cuss their interpretation. Finally, Sec. M contains our 
conclusions. 



degrees of freedom for each Mn spin. The band electron 
system and the local moments are coupled by a ferromag- 
netic exchange interaction that favors parallel alignment 
of the local moment and band electron spins. The total 
Hamiltonian of the system is: 



Ti, = Tin + TLt + Tis 



'T_/SS I 'T_/Sd 1 njdd 1 nj /I \ 



The first four terms in the right-hand side of the equa- 
tion describe a conventional superlattice system with 
many weakly coupled quantum wells: 

(i) Tio is the Hamiltonian for independent electrons 
in N isolated quantum wells. Its energy spec- 
trum is purely single-particle-like, and the quasi- 
particle spectrum is that of an isolated quantum 
weh 2DEG Ej{k\\) = Ej + S,{k\\), where k\\ is the 
wave vector parallel to the MQW heterointerfaces, 
^(fc||) = h'^k?J2m*, m* is the effective mass, and j 
is the quantum well subband index. We will take 
fc|| to be a continuous index , disregarding Landau- 
level formation in the weak magnetic fields we will 
consider. 

(ii) Ti-T contains the tunneling amplitudes that cou- 
ple quasiparticles in different quantum wells. In 
weakly-coupled superlattices it is a good approxi- 
mation to treat this term by leadin g or der pertur- 
bation theory as we discuss in Sec. II A. 



(iii) Tiscatt contains the scattering terms within a quan- 
tum well that allow a non-equilibrium quasiparti- 
cle to relax its excess energy, but does not contain 
terms that permit the quasiparticle system to bring 
its spin-subsystems into equilibrium. (These terms 
are absorbed in Hsi) Because it is difficult to de- 
scribe these scattering processes accurately, or even 
to know what they are in particular systems, we will 
use a phenomenological relaxation time approxima- 
tion. The time scale associated with thesaflrocesses 
is typically rather short (rgcatt '^ 0.4 ps).lij 

(iv) Hfjjj is the electron-electron interaction in the con- 
duction band for which we will use a Hartree mean- 
field approximation. (See Sec. |IIC|). 



II. THEORETICAL MODEL 

Following sucessfulO early work on bulk systems by 
Kossut,E3 Bastardp^ and Gaj,Ej we account for the pres- 
ence of Mn ions in DMS quantum wells, by combining 
a phenomenological exchange model with a virtual crys- 
tal approximation and mean-field-theory. The lattice pa- 
rameters and the band Hamiltonian parameters of a II- 
VI heterostructure are assummed to change smoothly as 
Mn++ spins are introduced in the system and a, S — 5/2 
quantum spin is assumed to be added to the low energy 



The remaining terms in Eq. (Q) describe spm-related 
physics: 

(i) Wfjfj is the exchange interaction between s con- 
duction band electrons and Mn local moments, an 
interaction that turns out to be ferromagnetic in 
II- VI MQWs. When the mean-field and virtual 
crystal approximations are employed, the effect of 
this coupling is to make the subband energies spin- 
dependent in those quantum wells that contain Mn 



ions: Ei 



EJ- 



(ii) 7i[5it represents the antiferromagnetic super ex- 
change interaction between Mn spins on neighbor- 
ing lattice sites that has been found to be important 
in modeUing bulk DMS systems .□ Since our inten- 
tion here is to address the qualitative physics of 
field-domains in DMS MQW systems, we neglect 
'^int- We do expect, however, that these interac- 
tions will be important for detailed modelling of 
specific experimental systems. 

(iii) Hsf contains the microscopic processes that allow 
equilibrium to be established between spin subsys- 
tems within a quantum well. The fact that spin r^ 
laxation can be quite slow in the conduction banded 
is one of the motivations for this work. Relaxation 
times in excess of 1 ns have been established ex- 
perimentallyS in II- VI semiconductor QWs without 
Mn. In II- VI DMS QWs these times are reduced 
to tens of picoseconds (but still larger than rscatt)-^! 
W e dis cuss the role of these terms at greater length 
inlilB. 



A. 2D-2D Tunneling 

The standard theory of tunneling relates the electric 
current between weakly-coupled subsystems to tunncliua 
matrix elements and subsystem spectral functions.Ej'EZI 
In our case we will apply this theory to describe the cur- 
rent flowing between one quantum well and its neighbor. 
Since elastic and inelastic scattering times in the quan- 
tum wells are shorter than any other time scale of the 
problem, we can follow the standard lines of tunneling 
theory and assume that the electrons in each well are in 
quasi equilibrium between succesive tunneling events and 
that their temperature is that of the lattice. We ignore 
interwell spin-flip processes, so that currents are carried 
between wells by the two spin subsystems in parallel. Ac- 
cordingly, the current per spin from the zth well to the 
{i -f l)st well is given by the following general expression: 









eV,) 



where a = 
dex, vq = 



X [/(^ - Mf ) - fie - A^f+i + eV,)] , (2) 

(t, J.) is the conduction electron spin in- 
is the 2D density of states per spin. 



kiki. 



2-Kn 



is the transmission coefficient between partic- 
ular wavevector states in the two quantum wells, eVi 
is the voltage drop across the iih barrier and /(x) = 
l/[eiq){xlkBT) + 1] is the Fermi factor, /i^ denotes the 
chemical potential in well i measured from the bottom 
of well i. A commonly used Lorentzian-shape function is 
chosen to represent the influence of disorder on quasipar- 
ticles in the jth subband within the zth quantum well: 



^S')-^ 



7 



^{e~£-{h)) +7^ 



(3) 



This form for the spectral function results from neglect- 
ing the real part of the disorder self energy, which intro- 
duces an unimportant rigid shift of the quasiparticle ener- 
gies, and the energy dependence of its imaginary part. 7 
is treated as a phenomenological parameter whose value 
may vary substantially from sample to sample and is to 
be taken from experiment. (7 — ?i/2Tscatt ~ 1 meV). In 
weakly-coupled superlattices the broadening due to scat- 
tering is much larger than the miniband width so that 
tunneling between quantum wells is sequential rather 
than band like. In this regime an electron undergoes 
many scattering events in one well before tunneling to 
the next well. Because of the epitaxial nature of the 
samples in question, we assume that the tunneling pro- 
cess conserves parallel momentum, i.e., effects such as 
interface roughness are not taken into account. This ap- 
proximation is made for the sake of definiteness and does 
not influence the qualitative conclusions we will reach. 
In addition, the typical electronic densities in particular 
quantum wells are assumed to be ^ lO^^cm"^ so that 
only the lowest subband is appreciably populated in the 
quasi-equilibrium state of a quantum well. Therefore, 
when current flows dominantly by a transition from the 
ground state of well i to the first excited state in well 
i -|- 1, rapid relaxation to the lowest subband (via, e.g., 
emission of a LO phonon) is assumed. This simplifica- 
tions lead us to an analytical expresion for Eq. (0) at 
r = OK: 






2n^h 



E^. 



27 



(i?fi-£;f+i^.+ey,)^ + (27)2 



(4) 



The sum in Eq. (||) is extended over all subands in 
the (i -I- l)st well. The transmission coefficients from the 
ground subband to the jth subband of the neighbouring 
well, Tj, are calculated by means of the transfer Hamilto- 
nian method.E£l Since this approach involves only orbital 
degrees of freedom, Tj does not depend explicitly on a. 
We are therefore implicitly restricting ourselves to the 
case when the spin splittings are not so large as to mix 
different subbands. The function S expresses the width 
of the energy "window" available for tunneling. For tun- 
neling between the lowest subbands: 



K - Ef, 


Mf+i - eV, < Ef, 


«+i-i?mi) 


/if + eF, < i?f+„ 


Mf - Mf+i + eV, 


otherwise , 



where Ei 1 is the first subband energy in the ith-well and 
so on. The expressions for the energy window for tunnel- 
ing to higher subbands are similar and differ only through 
the absence of Pauli blocking effects in the target layer. 



B. Spin splitting and spin relaxation 

The exchange interaction between the conduction band 
electrons and the Mn"*""*" ions produces a giant spin spht- 
ting of the 2DEG even in the presence of a small magnetic 
field. The bottom of the band for each spin rigidly shifts 
accordingly: 



EJ = E, - ,sA 



A . Js^NunSBs (^^) 



(5) 
(6) 



where Ej is the jth miniband center, s = +(— ) for 
a =T(i), Jsd is the exchange integral, A^Mn is the density 
of Mn^+ ions with spin equals to S* = 5/2, Bs is the 
Brillouin function, and TcS is an effective temperature 
which can, in general, include a correction due to antiferj-, 
romagnetic interactions between neighboring Mn ions.Q 
We have assumed here that the magnetic field (B) ori- 
entation is such that the quasiparticle energy is lowered 
for up spins. Note that our virtual crystal approximation 
for the Mn ions implies that the mean-field experienced 
by band electrons of spin a is spatially uniform. In this 
paper we take the field experienced by the local moments 
in Eq. (^ to be the external magnetic field. In fact, the. 
mean-field approximation we employ can be extendedES 
to include the contribution of spin-polarized band elec- 
trons to the total effective magnetic field experienced by 
the local moments. When this is done, ferromagnetism 
results. That extention of the models and approxima- 
tions employed here does in fact appear to account for the 
main ff^^l^ures of the carrier induced ferromagnetism that 



occursHJ LU in doped diluted ferromagnetic semiconduc- 



sJtion temperatures in these systems can 
^^' for p-type systems and for higher car- 



tors. The trar 
be substantial! 
rier densities. In the case of relatively low-density n-type 
systems, the ferromagnetic transition temperatures will 
be low and ferromagnetism does not necessarily occur, 
when direct interactions between the Mn local moments 
are included in the theory. Experimentally, most n-type 
DMS 2DEGs show no evidence for ferromagnetism. In 
this paper we will assume that the contribution of the 
spin-polarized conduction band system to the effective 
field experienced by the Mn local moments is negligible. 
If ferromagnetism did occur, the non-linear field-domain 
transport physics we discuss in this paper would be fur- 
ther enriched. 

In modelling spin-relaxation within a quantupi well we 



start by nelecting the transport currents Jj 



i,i+l 



For in- 
stantaneous spin-splittings smaller than the Fermi energy 
we use the following phenomenological spin-relaxation 
rate equation within each quantum well: 



dt 



f^i -f^i 



Tsf 



-1^0, 



(7) 



where a is the spin opposite to a and Tsf is the spin- 
scattering time. (We neglect the dependence of Tgf on 



A, which may be important in certain samplest3). This 
equation implicitly assumes the linear non-interacting 
2DEG relationship between density and chemical po- 
tential, |So that correlation effects are not taken into 
account £3 



M^^f-E^,). 



(8) 



In the absence of a driving bias voltage, the solution 
of Eq. (0) at i -^ 00 is nj = {ui + i^qA) /2 and n'l — 
{iii — t'oA) /2. Here rii — nl + ni is the total density of 
carriers in the ith well. We can see from this asymptote 
that Eq. (M) is valid only for spin splittings smaller than 
the chemical potential. 

For A greater than the chemical potential, Eq. (R) must 
be modified using Eq. (||): 



dt 
dt 



Tsf 

dt 



Mi 



E^ 



Tsf 



-I/O 



(9a) 



(9b) 



For large enough A, Eq. (||) leads to an equilibrium state 
with full spin polarization. (See Fig. H^). 




FIG. 1. Schematic illustration of spin-relaxation within 
a quantum well. The dashed lines indicate initial 

non-equilibrium chemical potentials while the dot-dashed line 
the equilibrium (t — > 00) chemical potentials. Panel (a) is for 
the case of spin-splittings smaller than fi^ — El while panel 
(b) is the case of for spin-splittings larger than fi^ — El for 
which the equilibrium state is completely spin polarized. The 
zero of energy in these plots is the electrostatic potential of 
the quantum well. 

Adding transport currents to these considerations 
leads to the following discrete continuity equations for 
the spin population in each quantum well: 



dnl 
dt 






Li+l 



P-i 



Mi 



Tsf 



-I/O 



!,■ 



,N 



(10) 



for the case /zj — E]^ > 2A. Otherwise, Eq. (^ must re- 
place the second term on the right-hand side of Eq. ([k]). 



C. Electrostatics 

In large area heterostructures, the Coulomb interac- 
tion is usually included in a Hartree mean-field approx- 
imation. The Poisson equation relates the electrostatic 
potential drop across MQW barriers, Vi, to the charge 
distribution among the quantum wells: 



V, - V,. 



{ui - N^) 



1,- 



.N. 



(11) 



where d is the superlattice period, e is the sample av- 
erage permittivity, and N^ denotes the doping density 
within the wells. (Experimentally, doping is usually ac- 
complished by placing a ZnCl2 layer in the barrier layers; 
this difference in electrostatics eopipared to our model 
has no important consequences. )liil 

By inspection of Eqs. (|I^) and ([ij) it is obvious that 
a set of boundary conditions must be provided for no 
and rijv+i- Within our model these layers play the role 
of source and drain, i:espectively. A simple way to repre- 
sent source and drainEj is to fix the density in both layers 
at high values. We take 



"■0, nN+i = KiVu 



(12) 



where k > 1 ^laji adjustable parameter. More sophisti- 
cated modelsP^rl have proven that a proper description 
of the contacts can have a strong effect on the selection 
of the the transport equation solution when multistabil- 
ity ocears, especially when dynamical solutions are al- 
lowedO However, we choose not to delve into these ef- 
fects in detail here since our main interest is on spin ef- 
fects. For the sake of definiteness, the contacts are taken 
to be unpolarized throughout our calculations; including 
spin-polarized injection in our theory would be straight- 
forward and indeed this may be a very interesting avenue 
to explore in future experimental and theoretical studies. 
With this representation of the source and drain, fixing 
the overall bias voltage. 



N 



^ = E^^' 



(13) 



i=0 



closes the set of equations. 



III. NUMERICAL CONSIDERATIONS 

Our model contains 5-/V + 1 unknown functions of time 
{2N chemical potentials, fil , 2N electronic densities, nf, 
and A'' -I- 1 voltage drops, Vi). These unknown func- 
tions are determined by the 2N constitutive equalities 
(Eq. d)), 2iV rate equations (Eq. (0)), N Poisson rela- 
tions (Eq. (p^)), and the total bias condition (Eq. (Oh). 
Thus, giving a physically sensible initial profile, the sys- 
tem of algebraic-differential equations can be integrated 
to yield a definite solution. Standard numerical methods 



are employed in solving Eq. (|l^) , being careful to use the 
appropriate spin-relaxation equation (either Eq. ([7|) or 
Eq. (^)) at each time step. 

The total current density J{t) traversing the sample at 
time t is determined by the following procedure. Differ- 
entiate Eq. (O) with respect to time and substitute the 
result into the sum of Eq. ( |lO| ) over spin indices in order 
to eliminate densities and chemical potentials. This leads 
to the following current which has the same value when 
evaluated for any well index, i: 



Jit) 



d dt 



+ Ju^+l{t) 



(14) 



The first term of the right-hand member of the previ- 
ous equation is the displacement current whereas the sec- 
ond term is the tunneling current, Ji^i+i{t) = JJ i^i{t) + 

J'ii+ii'^)- Foi' static steady state solutions discussed in 
this paper only the latter term is finite. 



IV. RESULTS 



We focus on ZnSe/Zni_^_yCd^Mnj,Se DMS MQWs. 
A value of x ^ 0.2 has been chosen to be consistent with 
barriers (~ 200 meV) sufficiently high to capture more 
than one subband in the quantum wells. For definite- 
ness, we focus on the case where Mn has been incor- 
porated only in the central well of the superlattice. In 
experimental samples the value of y can be varied over a 
wide range. We expect that the field-tuned field-domain 
effects we discuss will be strongest at moderate Mn den- 
sities, large enough to give rise to sizable spin-splittings 
but not so large as to increase the spin scattering rate 
excessively. 

The experimental samples reported on in Ref. Hq pos- 
sess ZnSe barriers too thick (^ 30 nm) for perpedicular 
transport. Under those conditions, the coefficients Tj 
are so small that electron tunneling would likely occur 
via impurity channels in the barriers or through different 
symmetry points in the host semiconductor band struc- 
ture, violating the assumptions of our theory. For these 
model calculations we choose a smaller barrier width: 
6 = 5 nm. The remaining phenomenological parameters 
have been fixed on the basis of available experimental 
data: it; = 10 nm, to* — O.IGtoe, A^^, = 2 x 10^^ cm~^, 
Tgf = 10 ps, and k = 1.5. w is the well width and k 
specifies the source and drain densities (sec Eq. (112)). 

In Sec. p^VA we study self-consistent steady-state solu- 
tions for low voltages. This will help us understand the 
more complicated behaviors that resul t from instabilities 
in the NDC voltage regime. (See Sec. IVB below.) 



A. Linear regime 

The behavior of the spin-dependent current density in 
Eq. (|j) depends on many variables. In Figs. || and 



for the case where 



Mi 



we have plotted Jf^+i, tor tne case wnere /ij, 
have been set to their equilibrium values (/ijjj=/i^ =/!,„). 
Thus the spin-dependent chemical potential, measured 
from the bottom of the well, will depend on the value of 
A (see Fig. |I|). 
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FIG. 2. (a) Current density flowing from a non-magnetic 
well to a magnetic one for A=0-6 meV in steps of 1 meV. 
The rightmost curve corresponds to A=0 meV. (b) Same as 
(a) but the carriers are now flowing from the magnetic well to 
the non-magnetic one. The leftmost curve is for A=0 meV. 

Fig. 0(a) plots the up-spin current from a non- 
magnetic well to a magnetic one, Jj„, i.e., net electron 
flow from a magnetic well to an non-magnetic one. For 
A = meV and low voltages the behavior is ohmic, as 
expected. A larger bias results in the alignment of the 
ground states of both wells within 7, giving rise to a first 
maximum in the current. After this bias is exceeded, 
Eq. (H) implies the appearance of a NDC region due to 
subband mismatch and of a second peak at higher biases 
when the first subband is aligned with the second sub- 
band of the next well. Increasing A decreases the value 
of the first peak since the up-spin density in the magnetic 
well increases and fewer states are available for tunneling. 
The peak corresponding to tunneling from first— >second 
subband is not affected in magnitude but its position is 
shifted to lower bias voltages because the bottom of the 
up-spin subband goes down. For a given value of the 
spin splitting (A '^ 3 meV), the magnetic well is fully 
polarized so that the first peak magnitude can no longer 
vary its value, a displacement due to subband lowering is 
observed instead. In Fig. 0(b) the up-spin current from 
a magnetic well to a non-magnetic one is shown. Here 
the first peak becomes larger as /xj„ grows since more 
electrons take part in the tunneling. This increase ceases 
once the fully polarized situation is achieved. In addi- 
tion larger voltages are needed to align the energy levels 
as A increases and the peak moves to higher voltages. 
Field-domain physics is at heart controlled by the in- 



terplay of electrostatics and the layer-to-layer non-linear 
I-V relationships. These examples illustrate that the 
layer-to-layer non-linear I-V^s can be altered by A (and 
hence an external magnetic field) when one of the layers 
contains Mn local moments. It follows that the field- 
domain structure must be A and field dependent as we 
show below. 
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FIG. 3. (a) Same as Fig. Q(a) for down spins. The leftmost 
curve corresponds to A=0 meV. (b) Same as Fig. Hb) for 
down spins. The rightmost curve is for A=0 meV. 



The corresponding down-spin current is plotted in 
Fig. ||. The considerations explained above for the up 
spin case may be invoked to understand this figure. For 
A > 3 meV (see Fig. ||(b)) the down-spin current flowing 
from the magnetic to the non-magnetic well (Jm„) van- 
ishes for any value of the applied bias since in the fully 
polarized regime, no down-spin carriers are present in the 
magnetic well. 

A numerical calculation for iV=3 gives rise to the 
MQW current-voltage characteristics shown in Fig. ||. 
A has been set to 3 meV. At low voltages the behav- 
ior is linear, up to the peak which marks alignment to 
within ~ 7 of the first subbands of all wells. Then, af- 
ter entering the NDC region, an increase of the current 
occurs because an electric field domain has formed. The 
underlying mechanism can be understood in the follow- 
ing terms. As the voltage is increased, the system prefers 
to mantain regions of low and high electric fields where 
intra-subband and inter-subband resonant tunneling can 
occur, rather than maintaining a constant field. The high 
field regime forms in the downstream side of the sam- 
ple. In this way, the superlattice minimizes the total 
resistance and current increases. In order to build the 
inhomogcnous field, according to Poisson equation (|ll| ) 
an excess electron density, or a monopole, must accumu- 
late in the well at the boundary between low and high 
field regions. As the voltage is further increased, config- 
urations with particular monopole locations sucessively 
become unstable and the monopole moves by one quan- 
tum well to increase the width of the high field region. 
Monopole motion gives rise to an abrupt decrease of the 
current, which is followed by a gradual increase as the 



voltage increases further. These successive monopole po- 
sition jumps lead ultimately to the sawtooth-like I-V 
curve shape of Fig. ^. The inset in this Figure shows 
the corresponding A = mcV I-V curve for compar- 
ison. The most obvious change is the appearance of a 
new branch that develops close to the first— ^second sub- 
band transition (sec below). 
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FIG. 4. I~V characteristics for A = 3 meV. The first peak 

occurs near where the first— +first subband resonant condition 

occurs in a configuration with no field domains. The following 

three branches reflect steady states with monopoles in one of 

the three quantum wells of the system. The monopoles move 

upstream with increasing bias voltage, increasing the number 

of high field barriers. Note that the final peak in the I-V curve 

is split in the spin-dependent transport case. Results for the 

A = meV case are shown in the inset for comparison. The 

vertical lines that connect different branches of the I-V curve 

are guides to the eye. 



The somewhat complicated behavior of the bias volt- 
age dependence of the steady state spin polarization is 
depicted in Fig. g for iV = 3 and A = 3 meV. Notice 
that in Fig. |5|(a) a non-zero steady state spin polariza- 
tion is induced in the non-magnetic wells by the spin- 
dependent transport currents. We have assumed here 
that the spin-relaxation time constant Tgf is the same 
in all quantum wells (10 ps). In practice these times are 
likely to be considerably larger in the non-magnetic wells, 
and the induced spin-polarization in these wells will have 
an even larger importance. For low voltages, the inset of 
Fig. ||(b) shows a reduction of the magnetic well spin 
polarization. This effect can be understood by realiz- 
ing that the minority spin current from a non-magnetic 
well to a magnetic one, J^„, grows when the system is 
driven out of equilibrium by a small bias voltage. This 
occurs also in jl„ but at a smaller rate because of the 
spin subband displacement in the magnetic well (com- 
pare Figs. |^(a) and y|(a)). Therefore, at low voltages the 
magnetic well polarization must decrease. Correspond- 
ingly, the upstream-well polarization must increase (see 
the inset of Fig. ||(a)) since down spins leak out of it. 
The enlargement of the downstream-well spin polariza- 



tion has the same origin, (see Figs. ||(b) and ^(b)) 
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FIG. 5. (a) Spin polarization for the upstream (down- 
stream) quantum well is represented by full (dashed) lines. 
The inset highlights the low voltage regime, (b) Same as (a) 
but for the magnetic quantum well. 



B. NDC regime 



As explained above, once the bias voltage reaches the 
NDC regime, electric field domains form in the sample. 
As we now discuss, their formation strongly influences 
spin-polarizations in both magnetic and non-magnetic 
quantum wells, with discontinuities associated with ev- 
ery break in the I-V curve. The magnetic well polariza- 
tion, P, varies particulary strongly, especially when the 
monopole moves through the magnetic well, and becomes 
stronger as the spin-splitting is increased. This can be 
seen in Fig. ^, which describes the general behavior of P 
with voltage and spin splitting. 
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FIG. 6. 3D plot describing the dependence of the magnetic 
well spin polarization on voltage bias and spin-splitting A. In 
constructing this^gure, Tst has been increased to 10 ns (still a 
reasonable value)El in order to magnify the effects commented 
on the text. 



If A is further increased (A=6 meV), new branches ap- 
pear in the I-V curve (see Fig. 0). The extra branch that 
appeared close to the Ei — > E2 transition at smaller A 
is now fully developed. In addition, a new branch forms 
before the Ei — > Ei transition. 



200- rn 




^"^'o 0.05 0.1 0.15 0.2 

Voltage (V) 

FIG. 7. I~V characteristics for A = 6 meV. Note the app- 
parance of new branches in the steady- state I-V curve. The 
spin polarization in the magnetic quantum well at low volt- 
ages is plotted in the inset. 



To better understand the subtle interplays that control 
these features, we have studied the self-consistent steady- 
state well-dependent spin-polarizations at the particu- 
lar voltages marked in Fig. 0. The series of MQW 
electrostatic profiles illustrated in Fig. g are dominated 
by classical field-domain physics not qualitatively influ- 
enced by the spin-dependent nature of the transport. 
Fig. P(a) (Vi) is the highest voltage at which intrasub- 
band {Ei — > Ei) resonant tunneling can be maintained. 
The electric field drops almost linearly along the system. 
In Fig. g(b) (V2) the formation of a high electric field 
domain in the last barrier is clearly observed, favour- 
ing a resonant condition between the third well and the 
collector first-excited subband.E3 The second branch in 
the NDC region (V3) involves the generation of a larger 
high field domain (see Fig. ||(c)). The domain wall is 
now located in the magnetic well. A jump to the first 
well (Fig. ^(d))(V4) is accompanied by further expansion 
of the high field domain. In this situation all tunneling 
within the wells takes place between the ground and the 
excited states, followed by a rapid relaxation to the first 
subband. 
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The spin polarizations plotted in Fig. H(a) are not observable 
on the scale of these figures. 



Fig. ^ illustrates spin-populations near I-V features 
where the spin-dependent element introduced by the 
magnetic quantum well plays a qualitative role. Fig. ||(a) 
describes the position of energy levels for the voltage V5 
marked in Fig. 0. The subband energies in the non- 
magnetic wells {E'^ and E^) are quasidegenerate. Notice 
that resonant tunneling occurs between E^ ^ and E2 1 ■ 
Further increase of the voltage, however, results in a de- 
crease of the current since now i?|'/ is then off-resonance 
(see Fig. ||(b)). The current is then increased again since 
i?{ I starts to match E2 1 ■ This explains why P does not 
show the behavior observed for smaller A (A=3 meV) in 
the linear regime. (In contrast to the inset of Fig. |^(b), 
the inset of Fig. ^presents a fiat polarization at low volt- 
ages.) P starts to increase only when V is such that _E|j^ 
reaches within 7 of £^2 1 • The splitting of the branch at 
higher bias voltage where all transport occurs via inter- 
subband {Ei —> E2) resonant tunneling has a similar ex- 
planation. The Fig. H(c) is stabilized by alignment of El j^ 
with E22 levels. Subband mismatch at higher bias gives 
rise to a sharp reduction of the current, which increases 
later, as the voltage is increased and resonant tunneling 
between Sfj^ and E22 levels is achieved. Incidentally, 
the dips in P illustrated in Fig. g can be explained in 
this way. When the latter alignment occurs, a large fiow 
of down-spin carriers streams towards the magnetic well, 
causing a sharp decrease of the polarization. The re- 
maining features in Fig. p^ can be understood in similar 
terms. 
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FIG. 9. Same as Fig. | for Vg, Ve, W, and Vg ((a), (b), (c), 
and (d) respectively). 



the one corresponding to the Ei — > £'2 resonance, but we 
choose not to show them here since these involve tran- 
sitions to higher excited energy levels, close to the top 
of the barrier, where our model breaks down. However, 
for sufficently high barriers it would be natural to obtain 
such behavior. 

Fig. ro shows three different values of spin polarization 
in the magnetic well which can be obtained at a partic- 
ular bias voltage, depending on the sample history, up 
sweep from zero voltage, down sweep from a high volt- 
age, and up sweep to an intermediate voltage followed by 
down sweep. We emphasize that this kind of hysteretic 
phenomena between magnetic states is driven here by 
electric fields. 



V. CONCLUSION 
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FIG. 10. Multistability of distinct polarization steady 
states within the magnetic well of a A'' = 9 MQW system. 
The inset shows a blow-up of three different steady states 
reached at 1/ = 0.08 meV. State labeled 1 (3) is achieved by 
sweeping voltage up (down) from a high (low) initial bias. 
State labeled 2 is obtained by sweeping voltage up to 1/ = 0.1 
(marked with a cross) and then reversing the sweep direction. 

So far all results were calculated by sweeping voltages 
up. Steady state solutions in the NDC bias voltage re- 
gion are in general multistable. We can obtain different 
solutions at a given bias V|©ltage by evolving solutions 
following different historiesj2l for example by decreasing 
voltages from a high initial bias. For a given voltage 
different values of the current with different density and 
spin polarization profiles may be achieved. To amplify ef- 
fects, in studying this possibility, we have set Tgf — 10 ns 
in a iV = 9 superlattice. The change in P with voltage is 
now so greatly increased (see Fig. |lO| ) that even reversed 
polarization can be observed at somewhat larger values 
of Tsf in the Ei — > £'2 resonant tunneling regime (not 
shown here). This is a direct consequence of the emer- 
gence of dominant spin bottlenecks. lJ The alterations in 
polarization are more visible for voltages greater than 



We have introduced and studied a simple model for 
growth direction non-linear transport in multiple quan- 
tum well systems containing magnetically doped layers. 
Our analysis is based on a tunneling Hamiltonian ex- 
pression for the current between spin-polarized quantum 
wells and on a phenomenological expression for spin- 
relaxation within quantum wells. Numerical studies of 
this model show that it predicts rich behavior due to its 
non-linearity and due to the additional degrees of free- 
dom introduced by spin-dependent transport. Nonlinear- 
ity manifests itself in the formation of electric field do- 
mains when the differential conductivity between neigh- 
boring layers is negative. We find that the spin polariza- 
tion of electrons in magnetic wells can change substan- 
tially when the system jumps between different branches 
of the I-V curve. When current flows finite electron spin- 
polarization extends from the magnetic quantum wells 
to non- magnetic quantum wells. When large equilibrium 
spin splitting, the I-V curve is strongly affected by the 
appearance of extra branches, due to tunneling into and 
out of spin polarized well subbands and these effects be- 
come more and more prominent when the characteristic 
time for spin-relaxation is longer. 

The effects addressed in this paper could be investi- 
gated experimentally by studying transport properties 
and by studying the polarization dependence of inter- 
band optical absorption and photoluminescence in MQW 
systems containing magnetically doped layers. The pre- 
dicted sensitivity of transport properties to external mag- 
netic fields, suggests that these systems could potentially 
be useful for magnetic field sensors, most likely in geome- 
tries with a relatively small number of quantum wells. 
The sensitivity to external fields will be strongest at low 
temperatures where the the Mn ions are easily polarized 
to produce the maximum equilibrium spin-splitting field. 
To illustrate the sort of effects that we expect to occur, 
we have considered only relatively simple geometries with 
a single magnetic layer. Other effects will occur in larger 
MQW systems with particular geometries. In general 



there is considerable lattitude for designing the Mn den- 
sity distribution in the MQW system to reahze desired 
magnetoresistance effects that could be described with 
the type of model we have introduced here. 

Analogs of the magnetotransport effects we discuss will 
also occur in ferromagnetic multiple quantum well sys- 
tems, similar to the delta-doped lasiered (Ga,Mn)As sys- 
tems studied by Kawakami et a/.cll These systems are 
ferromagnetic and the carriers are holes rather than elec- 
trons, leading to strain-sensitive spin-orbit-coupling in- 
duced magnetic anisotropyE3c3 and coercitivities. These 
properties suggest a rich interplay between the hysteretic 
magnetoresistance effects common in thin film itiner- 
ant electron magneta^a and the hysteretic effects dis- 
cussed here, which have their roots in electric field do- 
main structures. c3 Had ferromagnetism been taken into 
account in our model calculations, by solving for A self- 
consistenly in Eq. (|6|), not only the electric field domain 
structure but also the magnetic state configuration would 
have been sensitive to the bias voltage history. Exploring 
these possibilities appears to be a promising avenue for 
future experimental and theoretical work. 
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